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I. INTRODUCTION 


A. SEPARATED FLOWS 

The separated flow about bluff bodies has been almost completely unyielding to 
both analysis and numerical simulation for a number of mathematical reasons and 
fundamental fluid dynamic phenomena. Separation gives rise to the formation of free 
shear layers which roll up into vortex rings or counter-rotating vortices. They, in turn, 
Interact with each other, with the counter-sign vorticity generated at the base of the 
cody, and with the motion of often unknown separation points. The wake becomes 
unsteady even for a steady ambient flow and the problem of the determination of the 
characteristics of the wake becomes coupled to the conditions prevailing upstream of 
the separation points. Evidently, viscosity modifies radically the inviscid flow, which, 
in this case, cannot serve even as a first approximation to the actual flow. The 
boundary laver equations are not applicable beyond the separation points and are, 
therefore, of limited use in bluff-body flow problems. 

The separated unsteady flow situations involving wake return, as in the case of a 
decelerating or oscillating body, are an order of magnitude more complex than those 
where the vortices continuously move away from the body. The net effect of the wake 
return is twofold. Firstly, the proximity of the large vortices dramatically affects the 
coundary laver, outer flow, pressure distribution, and the generation and survival rate 
of the new vorticity. Secondly, the vortices not only give rise to additional separation 
points and/or additional vortices, but also strongly affect the motion of the primary 
vortices. These effects are further compounded by the diffusion and decay of vortices 
and by the three-dimensional nature of the flow. 

The existing finite difference and finite element methods cannot vet treat the high 
Reynolds number flows with sufficient accuracy for a number of reasons. The finite 
difference schemes require a very fine grid, a turbulence model, and a very large 
computer memory. It seems that the modelling of the turbulent stresses in the wake, 
particularly in time-dependent flows will be the major source of difficulty in all future 
calculations. Whether or not it will ever be practical to apply the finite difference and 
finite element methods to high Reynolds number flows is unknown. The inherent 
difficulties are certainly significant enough to warrant exploring other solution 


methods. 


1] 


Certain separated time-dependent flows may be simulated through the use of the 
discrete vortex model (DVM) {see ¢.c5 Chonn 1973) Sarpkava 1975). (iceicemoocen 
layers which emanate from the sides of the body are represented by an assembly of 
discrete vortices. The strength of the elemental vortices are determined through the use 
of the Kutta condition. The use of a suitable convection scheme enables one to march 
in time and to calculate the evolution of the wake, the velocity and pressure 
distributions, and the lift and drag forces acting on the body. The work described 
herein deals with the application of the DVM to decelerating flow about a two- 


dimensional camber. 


B. FLOW ABOUT A CAMBERED PLATE AND PARACHUTE COLLAPSE 

The determination of the deployment sequence of an axisymmetric porous 
parachute and the unsteady aerodynamic loads acting on it present a very complex 
coupled problem. The development of an analytical or numerical model which takes 
into account the effects of porosity, gaps, and variable opening schemes would allow 
numerical experiments on a large class of parachutes, reduce the number of expensive 
field tests to a few judiciously selected ones, and enable the designer to calculate the 
time history of the fall of the parachute and the strength required to survive the 
aerodynamic loads. However, the development of such a model is hampered by a 
number of difficulties. 

The previous models for parachute loads are based by and large on empirical 
assumptions (see e.g., Heinrich and Saari 1987; Mcwey 1972). They rely on the 
observation that families of parachutes open in a characteristic length and seem to 
have aerodynamic properties that relate well to the projected area of the parachute. 
The apparent mass 1s assumed to be a function of the projected area only and is not a 
function of the prevailing flow characteristics. The vortex sheet analysis was used by 
Klimas (1977) to derive the acceleration-independent apparent mass coefficient for 
arbitrary-shaped axisymmetric surfaces. Muramoto and Garrard (1984) used a 
continuous-source model to predict the steady-state drag of ribbon parachutes. The 
analvses did not, however, deal with the evolution of the unsteady wake and its 
interaction With the canopv. 

Itis in view of the foregoing that a fundamental study of the separated time- 
dependent flow about two-dimensional rigid cambered plates was undertaken. Clearly, 


the flow about a rigid cambered plate is considerably simpler than that about a porous, 


axisymmetric, and flexible parachute and the results, regardless of the degree of their 
agreement with corresponding experiments, may not have direct relevance to the 
practical problem under consideration. But the object of this investigation was the 
understanding of the evolution of the wake under controlled conditions rather than to 
provide a design tool. It is hoped that an investigation of this type will reveal the 
underlving physics of the phenomenon (particularly that of the parachute collapse), 
help to interpret the full-scale results and will provide inspiration for the development 
of more general vortex models with which the dynamics of axisymmetric, porous, and 


flexible parachute canopies can be investigated. 


Hi. ANALYSIS 


A. TRANSFORMATIONS AND THE COMPLEX YHVEGGinyY POIENT Es 

The calculation of the velocity of any one of the vortices and the force acting on 
the dDody requires a conformal transformation (in which the camber becomes a circle), 
a complex-velocity potential representing the vortices, their images, and the two- 
dimensional trrotational flow around the body, and the use of the generalized Blasius 
(NeEOrem 

The flow in the circle plane may be transformed to that about a cambered plate 
through the use of two successive transformations, one from § plane to the ¢ plane 


and the other from the >° plane to the z plane. Miese ane civem by (seen ian 


1)" 


a b. 


ee 





Figure 2.1 Cuircle and physical planes 


I4 


tJ 


oy 


(i ag 


and: €o= 4m (2.1) 


oO ’ 


YO 


Combining the two, one has a direct transformation from the § plane to the z 


plane as 


2 





Z =O. = za 


It is easy to show that the camber in the z plane 1s a circular arc. 

The y-axis in the z plane passes through the tips of the camber. It is 
advantageous to locate the origin of the coordinate axes at the geometric center of the 
Cee wee momsineme emer Ol ihe circle part Of whicly represents the camber. This is 


easily accomplished by shifting the origin of the coordinate axes by 


Mare Si a 
La is (225) 





Where Z,” is the x coordinate of the origin of the circle in the z plane. Thus, one has 


45 


— 





Z =e i 


v4 z = — >) 
Za cameo CAI 2 Z (2.4) 


Oo 


which transforms the circle in Fig. 2.la to the physical plane in Fig. 2.1d. Table 1 
summarizes the relationship between m, Zs the included angle of the camber, b, and 


the radius of the camber. 


s8 LES 
SUvIMIARY OF THE PARAMETRIC RELATIONSHIP 


m cies 2a b R= 1/m 
cos60 = .5 =| 120 866 2 
cos45 = .707 0. 180 .707 wld 
€0530 = 866 1/3 240 5 21/3 
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The complex potential function W im the Cirelemiaiaine (Scouse 
describes a uniform flow U (assumed to be time-dependent) with a doublet at the 
Origin to simulate the cylinder, ee clockwise-rotating vortices (called q-vortices), ee 
counter-clockwise rotating vortices (called p-vortices), and the images of all the p-and 


q-vortices in the circle plane may be written as 


PLn(S - =) 


ne a iP g 
We — CG — En (an ae 
é 2m P 2m Op 


m les i if 
te aelenG Sle) = ¥ Lk tng — =) —! Sy ENG — cgay 
eee k= 1°" Cp oe 
i, c? ee iP, i, : 
+ —M4in(g - =) - y ka (6 SG) a x q*Ln(g - (2.5) 
oi C =. q on 
“0g k= k oo 


in Which ae and Ck represent respectively the strength and location of the k-th p- 
VOTES re and Oka the strength and location of the k-th q-vortex, and c the radius of 
the cylinder; an overbar indicates a complex conjugate. The need for the separate 
identification of the p-and q-vortices and for the singling out of one of the vortices in 
Caclieshcar laver (iamely Pop and Vo, ie., the nascent vortices ) will become appurcnt 


later 


B. COMPLEX VELOCITIES OF VORTICES 

The convection of the vortices and the calculation of the forces acting on the 
body require the evaluation of the velocities at the vortex centers. Tor the velocities in 
the circle plane this reduces to subtracting from Eq. (2.5) the complexsiotencal 
corresponding to the vortex for which the velocity components are to be determined 
and evaluating the derivative of the remaining terms at § = ¢. To determine the 
velocities in the phvsical plane, however, one has to subtract (1 iV / 2) Ln(Z 92.) teem 


Eq. (2.5) or, in terms of € , the terms (Seeie cos any le crn 7) 


4 


K ) + Ln(l - oh (2.6) 
2m - a aa) 
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FoesioUlcmecmnetcd inal tiemiisa term in Eq. (2.6) is the complex function 


corresponding to the k-th vortex in the ¢ plane. The second term appears merely as a 


consequence of the transformation used. It is easy to show that for a p-vortex Eq. 
(2.6) may be reduced to (see e.g., Mostafa 1987) 


Ty, (7b, + m) 
Ae (Ga SEM SS aie 


(a7) 
A similar expression can be written for the q-vortex. 
©, 


KUTTA CONDITION 


The fact that the flow separates tangentially with a finite velocity at the edges of 
the plate (Kutta condition) may be expressed by requiring 





























dW 7 . 
— =0 a GS ia) eae (2.8) 
dy 
Thus, inserting Eq. (2.5) in Eq. (2.8) one has 
' iT, 1 l = Pog Pe l 
oT G. Sop _ Be — S0q a 
t + t 
SOp S0q 
m i. 
. ie ae I — spol bo 
Ka 27 WY, ¢ ee ee 4 faa 
a t é 
kp Kq 
% 
or 
el ees ent 


(2.9) 
Equation (2.9) may be decomposed into two parts as 


i 


to eegeee LV +) ame (2.10) 


where the terms containing the strength of the nascent vortices represent the velocity 
induced at the tip of the camber by the nascent vortices and the term in parenthesis the 
velocity at the tip due to all other vortices (and their images), the doublet at the center 
of the circle in the ¢ plane and the ambient velocity. 

Equation (2.10) represents two coupled equations for the strengths and positions 
of the nascent vortices. Thus, the solution of the said quantities does, in general, 
require an iteration. However, this iteration may be avoided by noting that the 
velocity induced by a nascent vortex at the opposite tip is very small and certainly 


negligible. Thus, Eq. (2.10) for one of the nascent vortices may be reduced to 


ee +(-u tiv.) = 0 (2.11) 
ae ig 2! 1V = re 
Pa Coa ok "a 
= 
Cog 


A similar expression may be written for the other nascent vortex. The use of the Kutta 
condition, as expressed by Eq. (2.11), will be further explained following the discussion 
of the tip velocity. It suffices to note that all nascent vortices satisfying the Kutta 
condition do not vield either the same tip velocity or the same velocity distribution in 
the neighborhood of the tip. There are, in fact, certain preferred positions for the 
nascent vortices which yield physically realistic velocity distributions near the tips of 
the cambered plate. These nascent vortex positions have been determined by Mostafa 
(1987). 


D. TIP VELOCITY 
According to the Kutta condition the tangential velocity at the tip is finite. It 
may be determined through the use of l’Hopital’s rule. 


The velocity at the tip is given by 


dw dW ad 
—=— — at z, = 
dz. dGwendz : 


| 
N 


_ = ib (2.12) 
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Vhesseonliaepital snulesthen yields: 


dw] d°W_ ib ae 
= ae =e (2.13) 
Z=Z, 


Equation (2.13) vields the desired finite tip velocities. 


FE. TIME DEPENDENT-FORCES 

The force acting on the body in the physical plane may be calculated either 
through the use of the pressure distribution or through the use of the rate of change of 
impulse. 


Bernoulli's equation for unsteady flow is given by 


ys 


2 Z 
a tL anes _ yh + >) -js = f(t) (2-74) 
p 2 p 7 





Where the indices indicate two points on the body in the physical plane. Since there is 
no pressure drop across the shear laver and since the integral term in Eq. (2.14) is zero 
at the tip (1.e., ds=0), one has 





oS Slt (2S) 


Where V,, and V,, represent the tangental velocities on the upstream and downstream 
faces of the tip. It is important to note that f(t) in Eq. (2.15) is also the time rate of 
change of circulation, 1.e., the rate at which vorticity is shed into the wake from the tip 
of the cambered plate. 

The normalized form of Bernoulli's equation between any two points m and n 


then becomes 


ie 


The integration of the differential pressure between the upstream and downstream faces 
of the camber yields the force components in the x and y directions, 1.e., the drag and 
lift forces. 

The force acting on the body can also be calculated through the rate of change of 


impulse. It is given by 


2 


ME, — 2.17 
) a a (Zee Ze (Qe) 


F = 4nxpc7U (1 — 
pc“U ( ae 





which may be written as 











ie ee 4 cU C _ m? 
a ~ 2pU2b_ a Uz yh oe? 
Cee Y @ e 
en ee 2.18 
ee ee ee (2.18) 


in which U, is the reference velocity; U, the rate of deceleration of flow and z= {(G,), 
i.e., the transformation given by Eq. (2.4). Equation (2.18) may also be deduced 
directly from the generalized Blasius equation. It is important to note that the force 
calculated from Ea. (2.18) ineludes the eflectwal the tate of chanee ol cicularen 
between two successive time steps. Thus, it may be smaller or larger (depending on the 
sign of 1) than the force calculated through the integration of the instantaneous 
differential pressure Eq. (2.16). This is because of the fact that the instantaneous 
pressure depends only on the prevailing flow conditions and does not account for the 
rate of change of total circulation between successive time steps. In the calculations to 


follow U, and c are taken as unity for sake of simplicity. 


F. METHOD OF CALCULATION 

The methods used in the past in the determination of the vorticitv flux from 
sharp-edged bodies may be roughly classified into two broad categories. The first of 
these involves the use of variable nascent vortex positions (see e.g., Sarpkaya 1968, 


1975) and the second, the use of fixed nascent vortex positions (see e.g., Clements 
1973-1975). 


In the present analysis the method of variable nascent vortex positions is used. 
To explain the computational details of the method let us consider a particular time t 
after the start of the motion and assumie t to be sufficiently large so that there are a 
number of vortices in the wake. Then the appearance and convection of the vortices 
proceed as follows: 


(1) Select a vortex position along the radial line defined by 8 = 117.72° (see 
Mostafa 1987). The very first location is taken r = 1.1. 


(2) Calculate the strength of the nascent vortex which satisfies the Kutta 
condition. This is an exact solution and requires no iteration. 


(3) Place the nascent vortex at the corresponding points in the circle and physical 
planes and calculate the tip velocity. 


(4) Calculate a new nascent vortex strength from ORSa - V,7)At where V, is the 
tip velocity and V,, the average of three velocities along the radial line in the 
Diysical plane, 1-e:,atr — 1-05,1,10 and 1.15. 


(5) Compare the newly calculated circulation with that obtained from the Kutta 
condition. If the difference between the two circulations is less than 0.001 
procede to the next step. If the said difference is larger than 0.001, carry out 
an iteration on the radial location of the nascent vortex as many times as 
necessary until the above condition is satisfied. If the circulation calculated 
from the Kutta condition is larger than that calculated from the tip velocities, 
the vortex must be moved towards the cylinder and vice versa. Also, each 
time the direction of the motion of the nascent vortex is changed (inward or 
outward), the marching distance is halved in order to accelerate the 
convergence of the two circulations. 


(6) Calculate the velocity induced at the center of all other vortices; 


(7) Convect the two nascent vortices with a velocity 0.5(V, + V, ) for a time 
interval At (note that the vorticity is convected with the average velocity of 
the shear layer). 


(S$) Convect all other vortices for the same time-interval At using a second order 
scheme given by 


2(t+ At) = z(t) + 0.5 [32(t) — 2(t—At)] At (2.19) 


in Which z = u + iv. 
(9) Remove the vortices from the calculation whenever they come nearer than 
0.05 to the camber in the physical plane (except the first 20 vortices from: the 


tip): 
(10) Coalesce the same sign vortices with a separation of less than 0.05 (in the 
physical plane, except the first 20 vortices); 
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(ll) Calculate the tangential velocities and pressures on the inner and outer faces 
of the camber. Determine the drag and lift forces through the integration of 
pressure and through the use of the rate of change of impulse. Make plots of 
suitable variables (e.g., velocity distribution near the tip, variation of nascent 
vortex circulation with time, evolution of the wake, etc.); 


(12) Check the flow conditions to determine the state of the calculations: 
(a) IfV,-V, > 0.1 repeat the foregoing steps; 


(6) Stop the introduction of nascent vortices if 0 < V, - V, < 0.1 and 
retumn to step UNone. 


(c) IfV, > V, switch the angular positions of the nascent vortices to their 
image points. Calculate V, as the average of the three velocities, at the 
upstream side of the tip of the camber, at three radial locations (0.95, 
0.9, and 0.85) and repeat the foregoing steps; and 


(13) Make plots of the variations of various flow parameters (e.g., tip velocity, 
nascent vortex circulation, evolution of the wake, force coefficients, etc.) and 
terminate the run. 


The foregoing steps are quite general and can be used for any camber, provided 
that the optimum points of placement of the nascent vortices are determined through 


the use of a method developed by Mostafa (1987). 
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Ili. DISCUSSION OF RESULTS 


A. NUMERICAL AND PHYSICAL EXPERIMENTS 


The calculations were carried out for a time-dependent normalized velocity given 





by 
U Cy | 
—-=1 for T= < 9.72 (3.1) 
ioe c 

and 
u a, 7 
— = 0.97T -—0.05T* — 3.70 (32) 
ue 

and 
A = 0.97 -0.10T (3.3) 


in the interval 9.72 <T <11.48 and 


L = —0.3423T + 0.0072T? + 3.82 (3.4) 
and 
A = —0.3423 + 0.01445T (3.5) 


Pi eemuesvainii.cs = |.=1/,.95. for I larger than 16, the velocity and acceleration 
are zero. These velocities and accelerations correspond to that encountered in a series 
of experiments carried out in a vertical water tunnel. A detailed description of the 
equipment and procedures is given by Sarpkaya and Ihrig (1986) and will not be 
repeated here. Evidently, the calculations can be carried out for any specified variation 


of the velocity. For the case under consideration. the flow begins to decelerate at T* 


ty 
we) 


= 9.72 and the velocity of the ambient flow reduces to zero at about T* = 18, (see 
Pie, 3.1) 

The computer program provided, at times specified, the positions of all the 
vortices, the rate of shedding of vorticity from the tips of the camber, the velocity 
distribution on the upstream and downstream faces of the camber, the differential 
pressure distribution. and the force coefficients. 

Figures 3.2 through 3.5 show, at T* = U,t/e = 6.00, the evolution of the wake, 
the differential pressure distribution, and the velocities at the upstream and 
downstream faces of the camber. Figures 3.6 through 3.17 show similar plots of the 
wake, pressure, and velocity at larger times. These and other figures (not reproduced 
here for sake of brevity) show that the characteristics of the flow develop symmetrically 
prior to the onset of deceleration (1* < 9.72) and the differential pressure is positive 
everywhere (1.e., the pressure inside the camber is larger than that outside). 

Following the onset of deceleration (see ¢.g.. Figs, 3.8, 3:12, .andegeic; une 
differential pressure near the axis of the camber becomes increasingly negative. The 
significance of this result 1s that had the model been flexible (as in the case of a 
parachute) the central part of the camber would have collapsed as a result of the 
particular deceleration it 1s subjected to. 

For IT* larger than about 13 (for the ambient flow under consideration), the 
velocities induced at the downstream edges of the camber by the large vortices moving 
sideways and towards the camber give rise to oppositely-signed vorticity. This, in turn, 
leads to the rapid growth of the secondary vortices (see e.g., Figs. 3.14 and 3.15 at T* 
= 21.275). Attention is drawn to the fact that the differential pressure shown in Fig. 
3.16 1s negative over a large central portion of the camber. 

Figure 3.18 shows the variation of the circulation of the nascent vortex with T*. 
As expected, the vorticity flux is quite large at the start of the motion. It decreases 
gradually prior to the onset of deceleration and rapidly thereafter. For T* larger than 
about 12, 1.e., when the secondary vortices come into existence, the circulation becomes 
negative. It reaches a minimum value when the ambient velocity and acceleration 
reduce to zero (at about T* = 18). The subsequent motion of the primary and 
secondary vortices increases the tip velocity and hence the strength of the nascent 
vortices. Nevertheless, the circulation remains negative, i.e., only the secondary 


vortices continue to receive additional vorticity. 


24 


Figure 3.19 shows the velocities V, and V, as a function of T*. The tip velocity 
V, decreases from an initially large value of about 3 to a nearly constant value of 
about 1.4 just prior to the onset of deceleration. Subsequently, V, decreases rapidly 
during the period of deceleration and prior to the inception of the secondarv 
separation. [Then V, increases to about 2 because of the backward motion of the large 
vortices near the tips of the camber. Finally, V, decreases once again as the primary 
and secondary vortices move sidewavs and awav from the tips of the camber due to 
their mutual induction (see Fig. 3.15). 

The variation of V, with T* is significant only during two, relatively short, time 
intervals: at the start of the motion and at the start of the deceleration. These are the 
periods during which the vorticity flux changes rapidly in order to maintain the Kutta 
condition. During the remainder of time V, is negligibly small, as expected on the 
basis of the pioneering experiments of Fage and Johansen (1928) with steady flow over 
various types of bluff bodies. 

Figure 3.20 shows the variation of the drag coefficient as a function of T*. It is 
calculated through the the integration of pressure about the camber. Figure 3.20 also 
shows that C, rises rapidly (due to the rapid accumulation of vorticity in the growing 
vortices) and begins to decrease as the vortices develop under the influence of a 
constant ambient velocity. Then the force decreases sharply at the onset of 
deceleration and goes through zero near the middle of the deceleration period (1* = 
14). The force acquires its largest negative value towards the end of the deceleration 
period. Subsequently, the force gradually decreases to Zero. 

Figure 3.21 shows a comparison of the calculated and measured drag coefficients. 
In general the agreement between the calculated and the measured drag coefficient ts 
quite good. The differences are primarily due to the fact that the diffusion of vorticity 
nas not been taken into account in the numerical analysis. I[t is possible to introduce a 
smiall artificial reduction in circulation in order to bring the calculated and measured 
values into closer agreement. This has been avoided in the present analysis 1n order to 
keep the discrete vortex analysis as pure and simple as possible. 

Finally, Fig. 3.22 shows a comparison of the normalized drag forces for various 
periods of the initial steady flow, prior to the onset of deceleration. It 1s clear that the 
drag forces bevond the point of deceleration are nearly identical. In other words, the 
force acting on the camber is not materially affected by the duration of the ambient 


steady flow within the range of the parameters encountered in the present study. 
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Figure 3.1 Variations of the velocity and acceleration with T* 
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Figure 3.2 Position of vortices at T* = 6.00 
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Figure 3.3. Instantaneous velocity of the vortices at T* = 6.00 
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Figure 3.4 Differential pressure distribution at T* = 6.00 
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Figure 3.5 Velocity distribution on both faces of the camber at T* = 6.00 
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Figure 3.7 Instantaneous velocity of the vortices at T* 
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Figure 3.8 Differential pressure distribution at T* = 12.275 
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Figure 3.9 Velocity distribution on both faces of the camber at T* = 12.275 


34 


Q°sg 


* ee 
oe € oi, ise 
e «? . 
e re * eee! e 
e 
°% ° i ie. 
e « 
e e° "ee é 
es e e S P) en 
e e * « 
® 
. e 
e ° et@ue - ° e e A 
e O.- ie ae e 
S ¢ e * ® e° oe be « e « e 
e © oe ° e re e 
e i 2 ee 
bid = ; e * 6 be S a. eet ae 
» e e e . a a A « - e 
ee ee? ° aa 
~ e : os e 
ie ° z ® ey fa : < 
e e e e = a e e 
e e e ¢@ ee A e 


We475 


Figure 3.10 Position of vortices at T* = 
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Figure 3.12 Differential pressure distribution at T* = 13.775 
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Figure 3.13 Velocity distribution on both faces of the camber at T* = 13.775 
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Figure 3.14 Position of vortices at T* = 21.275 
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Figure 3.15 Instantaneous velocity of the vortices at T* = 21.275 
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Figure 3.17 Velocity distribution on both faces of the camber at T* = 21.275 
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Figure 3.18 Variation of the circulation of the nascent vortex with T* 
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Figure 3.19 Variation of the velocities 


4 


0°S2 


S°c¢ 


0°02 


S41 


O°st 


J/ IN 
S*2l 


8 
ae OBE Gece 08002069998 0080008 08 208 0 0 08 beleterarer tere Sot CERARRSS CENLEENEN SENRERNEES CARO E ES BIDE E DEON ons oe va 
° 


O'Ol 


S°Z 


0°S 


s ¢ 





Q°O 


Q Calculated drag coefficient as a function of T* 


A 
~ 


onine 3. 


45 


Gz O's, 


Figure 3.21 


a” 
20.0 22.5 25.0 


725 


PHM 6 
15.0 


oe pr 
73 10.0 2.9 
UT/C 


el) 


eno 


9 
° 
0.0 


Q°S Gr Sas Oe, O°} 0°O Ot=— Uee= Ore 0 r- 


Comparison of measured and calculated drag coefficients 


46 


40.0 aa0 


ae 
=F 
joe0 


<= 
15.0 20.0 25.0 30.0 
UT/C 


10.0 


Sa0 


G0 


az a°9 O70 OTI- O°Z- O°E- O°F- 





Figure 3.22 Measured drag coefficients for various periods of initial steady flow 
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B. CONCLUDING REMARKS 

The results presented in this section have shown that the discrete vortex model 
can be used with confidence to predict the evolution of the wake about a cambered 
plate immersed in an arbitrary time-dependent flow. The drag coefficients resulting 
from the analysis and experiments agree reasonably well. This agreement can be 
improved with the introduction of a small circulation dissipation. 

The development of negative differential pressures near the central region of the 
camber 1s thought to be primarily responsible for the inception of the partial collapse 
of a parachute at high rates of deceleration. This phenomenon takes place even when 
the total drag force acting on the parachute 1s still positive. The sample analysis 
presented herein also shows that the negative differential pressure can cover a large 
region of the parachute and even result in negative drag. The basic idea emerging from 
the analysis reported herein is that the designs which incorporate into them the idea of 
delaying or preventing the return of the shed vortices to the canopy (e.g.,porosity 
management, change of deceleration history, parachute shape, dissipation and/or 
destruction of the organized wake) will be the ones which could avoid the collapse 
phenomenon. Extensive analysis and small scale experiments coupled with a few 
judiciously selected field tests may help to arrive at practically and phenomenologically 


sound parachute designs. 
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